Disorder-Induced Static Antiferromagnetism in Cuprate Superconductors 
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Using model calculations of a disordered d-wave superconductor with on-site Hubbard repulsion, 
we show how dopant disorder can stabilize novel states with antiferromagnetic order. We find that 
the critical strength of correlations or impurity potential necessary to create an ordered magnetic 
state in the presence of finite disorder is reduced compared to that required to create a single isolated 
magnetic droplet. This may explain why in cuprates like La2- :c Sr a; Cu04(LSCO) low-energy probes 
have identified a static magnetic component which persists well into the superconducting state, 
whereas in cleaner systems like YBa2Cu30e+«(YBCO) it is absent or minimal. 

PACS numbers: 74.72.-h,74.25. Jb,74.81.-g,74.25.Hc 



The occurrence of high-temperature superconductivity 
in cuprates near the antiferromagnetic (AF) phase of the 
parent compounds has prompted speculation since their 
discovery that superconductivity and magnetism were in- 
timately related. For the most part, it has been assumed 
that the two forms of order compete and do not coexist, 
consistent with the vanishing of the Neel temperature Jjv 
before the onset of a superconducting critical tempera- 
ture T c , and the suppression of T c near doping x = 1/8 
where static stripe phases can be stable On the other 
hand, there have been persistent reports of static AF at 
low temperatures T in the superconducting phase at low 
doping, as measured by muon spin resonance (/iSR) 0, Q , 
nuclear magnetic resonance (NMR) 0, @| , and elastic neu- 
tron scattering (NS) @, 0] . The NS experiments reveal an 
incommensurate (IC) ordering wavevector evident by a 
quartet of peaks surrounding (tt,it). Since the neutron 
response is enhanced by an applied magnetic field 0,0,0 
several authors have discussed it in terms of coexisting 
d-wave superconductivity (dSC) and field-induced spin 
density waves floj]. Recent magnetic Raman scattering 
data on LSCO has been discussed in terms of such ef- 
fects as well ll| . However, static order also exists at 



zero field in the underdoped phase of LSCOjl, Q, and 
has been attributed to disorder Q. In optimally doped 
LSCO, Kimura et aJ. jT3| did not detect ordered moments 
in pure and 1% Zn-substituted samples. An elastic peak 
similar to the pure underdoped material was observed 
when 1.7% Zn was added, however [l2j|. 

Phenomena similar to those in LSCO have been ob- 
served in other materials, e.g. Yi_. r Ca :E Ba2Cu306, and 
Bi 2 Sr 2 CaCu208+2; (BSCCO) where /J.SR directly reveals 
a slowing down and subsequent freezing of spin fluctua- 
tions as T is lowered 0, [l^. On the other hand, experi- 
ments on optimally doped YBCO, even with significant 
percentages of Zn, have never detected static magnetic 
signals. While there have been reports of AF coexist- 
ing with dSC in underdoped YBCO, recent NS measure- 
ments on YBCOg.s found that AF order, while static 



from the point of view of NS timescales, 10 s, was fluc- 
tuating faster than the timescale, 10~ 6 s, for //SR. Thus, 
it appears that while in YBCO low-frequency AF fluctu- 
ations are present, they do not "freeze out". Recently, 
reports of static magnetism in this system near O con- 
tent 6.35, close to the onset of superconductivity, have 
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been reported by ^SR[1J, |15J, but are still controversial; 
the main point is that the spin-glass phase is minimal in 
YBCO compared to LSCO[16[. There are many differ- 
ences between YBCO and the other cuprates, of course, 
but most importantly YBCO appears to be the cleanest 
material because the O dopants can order in the CuO 
chains. LSCO, on the other hand, is doped by randomly 
located charged Sr ions only 2.4 A away from the Cu02 
planes. Therefore, it seems likely that disorder itself may 
be responsible for inducing the magnetism in zero field 
in LSCO, and possibly part of the magnetic field depen- 
dence as well. 

/iSR experiments [3j indicate that several changes oc- 
cur with decreasing T for a given sample: one, where 
random freezing of moments occur, and a second, where 
a sharply peaked magnetization distribution arises. We 
propose that the first is due to the creation of isolated 
AF droplets, and the second to the formation of networks 
of such states which exhibit quasi-long range order. The 
latter is a subtle process, which, even for a single pair 
of impurities, depends on the relative orientation and 
distance of impurity positions 17 1 . However, the basic 
physics of ordering is easy to understand. In one dimen- 
sion Shender and Kivelson[18j pointed out that the inter- 
actions between impurities in a quantum spin chain are 
non- frustrating: if an impurity creates a local AF droplet, 
a second one can always orient itself to avoid losing ex- 
change energy. In two dimensions (2D) this continues to 
apply for spin models with nearest neighbor exchange, 
but may break down in the presence of mobile charges. 

Below we study the disorder-induced magnetization in 
a dSC with correlations described by the Hubbard model 
treated in an unrestricted Hartree-Fock approximation. 
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We focus on important qualitative differences between 
the clean and dirty limits of both the underdoped and 
optimally doped regimes. Specifically, we show that with 
a fixed choice of realistic parameters, as in LSCO, the 
magnetism is present at low doping, disappears at opti- 
mal doping, but can be recreated with a small concen- 
tration of strong scatterers. For YBCO such effects are 
absent since the disorder potential in the CuC>2 planes is 
negligible and independent of doping. This result agrees 
with a recently proposed origin of the unusual transport 
measurements in LSCO compared to YBCQ(l9ll20j. 
Model. The Hamiltonian, defined on a 2D lattice, is 

H = - 53 + \ Un i-" + v i mp - m) c\ a c ia 



E( A *A t T4i+ H - 



(i) 



(ij) 



Here, cJ CT creates an electron on site i with spin er, and 
tij = {t,t',t"} denote the three nearest neighbor hop- 
ping integrals, V- mp — Ylj=i V" np Sij is a nonmagnetic 
impurity potential resulting from a set of N point-like 
scatterers of strength V lmp , fj, is the chemical potential 
adjusted to fix the doping x, and Ajj is the d-wave pairing 
potential between sites i and j. The amplitude of A^ is 



set by the dSC coupling constant Vd|2lJ. We fix the band 
t' = — OAt and t" — 0.12i, giving the Fermi surface shown 
in Fig. DJa). We have solved Eq.(TT]) self-consistently by 
diagonalizing the associated Bogoliubov-de Gennes equa- 
tions in the T — limit on 34 x 34 systems 21]. 

Eq.JT]) has been used extensively to study bulk 
competing phases, field-induced magnetization, as well 
as novel bound states at interfaces between AF and 
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with considerable 
in Li- 



and Zn- 



superconductors [231] . It has also been used to study field 
induced moment formation around nonmagnetic impu 
rities in the correlated dSC 17 , 2^ 
success at fitting NMR lineshapes 
substituted YBCO. Interference effects of these param- 
agnetic states, while significant, are less important than 
in the heavily disordered case studied in this paper where 
the ground state in zero field exhibits local magnetism. 

Results. For clean systems it is well-known that there 
exists a large region in (U, Vd, T)-space dominated by spin 
and char ge o rdered stripe states of coexisting dSC and 
AF order [23j. For fixed Vd and T we denote by U C 2 the 
critical Coulomb repulsion for entering this bulk mag- 
netic state, i.e. for U < U C 2, in the absence of disorder, 
the ground state is a homogeneous dSC. 

For a single nonmagnetic impurity [3], there exists a 
lower critical U c % such that local impurity-induced mag- 



netization exists for U c \ < U < f/ C 2 [25| as shown in Fig. 
[TJb) . In general U c ± depends on the impurity strength 
V mp [l3,|23|- However, already for two impurities inter- 
ference effects modify the 1-impurity phase diagram [17 1. 
The magnetic structure factor S(q) associated with the 
AF droplet in Fig. [ljb) is shown in Fig. [ljc). It is 



dominated by four IC peaks along the diagonals, a result 
which is, however, sensitive to the input parameters. 

We now turn to the many-impurity situation. In the 
following model for LSCO the Sr ions are assumed to 
be the primary source of disorder, such that n lmp = x, 
where n mp denotes the impurity concentration. These 
systems are in the strongly disordered regime where the 
Fermi wavelength Xp is comparable to the average dis- 
tance between the dopants, such that the disorder is far 
from the 1-impurity limit and novel emergent spin-glass 
states can arise. Since the Sr dopants are removed from 
(but close to) the Cu02 planes we model them as inter- 
mediate strength scatterers with V lmp — 3.0t (below we 
use V dop to indicate the Sr potential). Fig. [2a ) shows 
the simplest schematic picture of disorder stabilization of 
a single Neel phase [18| by nonmagnetic impurities. Not 
all impurities in the correlated system need "magnetize" 
for a given U, however: in the disordered system, the 
effective criterion to drive the impurity through the lo- 
cal magnetic phase transition is different for each im- 
purity. Increasing the repulsion U then increases the 
concentration of impurities which induce a local mag- 
netization droplet, as shown in Fig. [2jb,c). With further 
increase of U, the system evolves from a state with dilute 
non-overlapping AF droplets to connected spin textures 
(Fig. EJd)). The resulting patterns exhibit AF domain 
structure and are more complex than suggested in Fig. 
[2ja), due to frustration induced by the charge degrees of 
freedom in Eq.JT]), and/or glassy supercooling effects. 

While the current mean field treatment of the Hubbard 
model docs not faithfully capture the band narrowing due 
to correlations near half-filling which leads to the Mott 
transition, it may be expected that underdoped systems 
are characterized by larger effective interactions. In our 
picture for LSCO, the x dependence of the spin order is 
therefore described qualitatively by the sequence [2fd) — > 
I2I C ) —* IHb), until it disappears completely at effective 
C/'s below U c i near optimal doping. Increasing x should 
also be accompanied by a weakening of the Sr potential 
ydop j ue £ enhanced screening. Within our model, in- 
creasing U or V dop leads to qualitatively similar results, 
and we cannot determine from this approach which ef- 
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FIG. 1: (Color online) (a) Fermi surface with x — 7.5% for the 
band used in this paper. For this band the critical value for 
bulk order is U C 2 = 3.50t. (b,c) Magnetization and structure 
factor S(q) for a single point-like scatterer V lmp — 3.0t for 
x = 7.5% and U = 3.3t. For these parameters, U c i = 3.25t. 
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FIG. 2: (Color online) (a) Schematic: stabilization of single 
Neel phase by impurities, (b-d) Disorder-induced magnetiza- 
tion for a single fixed impurity configuration with x = 7.5% 
and U = 2.8* (b), U = 3.2* (c), and U = 3.6* (d). In (a-d) 
the black dots indicate the impurity/dopant positions. 
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FIG. 3: (Color online) (a) Electronic density and (b) dSC 
gap for the parameters shown in Fig. [2jc). (c,d) display 
the impurity averaged magnetic structure factor S(q) for the 
dilute magnetic droplet limit (small U) and dense limit (large 
U) corresponding to Fig. [2jb) and Fig. EJd), respectively. 



feet is dominant in real systems. Note from Fig. [2] that 
AF droplets are induced for U > 2.4t, a substantially 
reduced critical value compared to the 1-impurity case 
in Fig. [TJ This is because the Hubbard correlations in- 
duce charge redistributions which alter the effective local 
chemical potential, such that the criterion for the mag- 
netization of each impurity depends on its local disorder 
environment. Some regions containing impurities have 
charge densities closer to the phase boundary for AF or- 
der, thus enhancing local moment formation relative to 
the single impurity case. In the limit of large U , the mag- 
netic order becomes qualitatively similar to that arising 
in a stripe state with quenched disorder [26j. 

In Fig. [3Ja,b) we show the inhomogeneous electronic 
density and dSC gap corresponding to Fig. |2^c). The re- 
pulsive Sr dopants locally suppress both quantities which 
are therefore anti-correlated with the induced spin order. 
As expected, the gap varies more smoothly compared to 
the density due to the coherence length £ of the dSC 
condensate. The magnetic structure factor S(q) (aver- 
aged over 30 distinct impurity configurations) associated 
with the disorder induced magnetic order in Fig. [2fb,d.) 
is shown in Fig. [3jc,d), respectively. It is dominated by 
an IC ring surrounding (n, 7r) with an intensity distribu- 
tion similar to the single impurity S(q) for the dilute case 
(Fig. [3jc)), but which rotates into a + shaped pattern 
for the connected spin textures at larger U (Fig. Etd)). 
This implies a rotation (and weakening) of the IC pattern 
close to the region where the static order disappears. As 
mentioned above, such details are, however, sensitive to 
e.g. the specific band parameters 27 1. 

Lastly we discuss optimally doped LSCO where static 
AF is absent in nominally clean samples, but where it 



can be induced by magnetic fields [6j or critical concen- 
trations of strong scatterers. For instance, Kimura et 

a] .[12| found that for x = 0.15 it takes approximately 
2% Zn to induce IC peaks in the NS diffraction. Below 
this critical concentration there was no measurable signal 
above the background. We have simulated this situation 
by solving Eq.((T]) with x — 0.15 in the presence of 1% and 
2% randomly distributed strong scatterers (V Zn = 100*) 
in addition to the n lm P = 15% weak (V dop = It) dopant 
impurities. As mentioned above the weaker V dop com- 
pared to the underdoped case is expected from an en- 
hanced screening of the Sr potential at optimal doping. 
For a single Zn impurity, magnetization is induced only 
for U > 3.35*. The many-impurity results including both 
ydop an( j yZn are snown f or different Coulomb repulsions 
in Fig. |4]and Fig. El For systems with 1% Zn Fig. [4|a- 

b) reveals a negligible induced magnetization. However, 
as seen from Fig. SJc-d), simply adding enough strong 
scatterers can induce sizable local magnetic order. It is 
remarkable that such a small difference in the concentra- 
tion can induce magnetism similar to the experimental 
observations. The number of Zn ions that induce AF 
droplets depends on the Hubbard U as seen by compar- 
ing Fig. 2] and Fig. O The resulting disorder-averaged 
S(q) agrees well with the NS measurements in Ref. \l% 

Conclusions. The interplay of dopant disorder 
and electronic correlations can induce novel magnetic 
states which we propose exist in intrinsically disordered 
cuprates like LSCO and BSCCO, in contrast to the 
cleaner YBCO system. Lastly, we showed that small 
concentrations of Zn can induce a similar magnetic state 
which has apparently been observed in NS experiments 
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FIG. 4: (Color online) Induced magnetization (left) and im- 
purity averaged S(q) (right) for U = 3.2t with 1% (top) and 
2% (bottom) concentration of randomly distributed strong 
scatterers. Note that the figures are shown on the same scale. 
For clarity, the black dots (left) only show the Zn positions. 
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FIG. 5: (Color online) Same as Fig. g]but for U = 3 At. 

on optimally doped LSCO. An obvious question is to 
what extent the magnetic state influences the scattering 
of quasiparticles in the dSC and normal states. Studies 
along these lines are in progress. 
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